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abstract. We consider the equation 



(1) e 2 Au = (u- a(x))(u 2 - 1) in Q, — = on dCl, 

where SI is a smooth and bounded domain in R", v the outer unit normal to dfl, and a a smooth function 
satisfying — 1 < a(x) < 1 in fl. We set K, £1+ and Q_ to be respectively the zero-level set of o, {a > 0} 
and {a < 0}. Assuming Va ^ on K and b ^ on 90, we show that there exists a sequence Sj — > 
such that equation JTJ) has a solution u ej which converges uniformly to ±1 on the compact sets of f2± 
as j — > +oo. This result settles in general dimension a conjecture posed in [19], proved in [15] only for 
n = 2. 
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1 Introduction 

Given a smooth bounded domain SI of R™ (n > 2), we consider the following problem 



, , J e 2 Au = h(x,u) in SI, 

X ' ( ' 11-0 on dQ, 



dv 

where e is a small parameter, the unit outer normal vector to <9S1 and h a smooth function such that 
the equation h{x 1 1) = admits two different stable solutions t\ 7^ t-z for any x £ SI. Using matched 
asymptotics, Fife and Greenlee in [TH] proved under some hypothesis on h the existence of a solution of 
© which converges uniformly to U in the compact subsets of Sl^, i = 1,2, where Qi and SI2 are two 
subdomains of SI such that SI = Slj U SI2. 



In this paper we consider a model heterogeneous case h{x, u) = (u — a{x))(u 2 — 1), for a smooth function 
a satisfying —1 < a(x) < 1 on SI and Va 7^ on the set K = {a(x) = 0}, with K n d^l = 0. We prove 
the existence of a new type of solution of © for any n > 2 settling in full generality a result previously 
proved in [15] for the particular case n = 2. 
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2 E-mail addresses: malchiod@sissa.it (A. Malchiodi), 
3 E-mail addresses: wei@math.cuhk.edu.hk (J. Wei) 
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Let us describe the result in more detail: in the case h(x, u) = (u — a{x)){u 2 — 1) problem ([2j> becomes 



(3) 



£ 2 Au=(«-a(j;))(!i 2 -l) infi, 

|2i = o on on. 



In particular, when a = 0, (J3J is nothing but the Allen-Cahn equation in material sciences (see [6]) 

. . J e 2 Au + u — u 3 = in fi, 

^ \ = on <9Q. 

Here the function u(x) represents a continuous realization of the phase present in a material confined to 
the region f2 at the point x. Of particular interest are the solutions which, except for a narrow region, 
take values close to +1 or —1. Such solutions are called transition layers, and have been studied by many 
authors, see for instance [IllEllIilClllHfflEDlia references therein 

for these and related issues. 

In this paper, we are interested in transition layers for the heterogeneous equation (J3J). Define 

K = {x e fl : a(x) = 0} . 

We assume that K is a smooth closed hypersurface of f2 which separates the domain into two disjoint 
components 

(5) o = o_ uKun+, 

with 

(6) a(x)<0inf2_, a(x) > in fi + , Va ^ on K. 
We then define the Euler functional J £ (u) associated to (J3|) in S7 as 

(7) Je{u) = £ - ( |Vw| 2 + f F(x,u)dx, 
where 

F{x 1 u):=j (s- a(x))(s 2 - l)ds. 

The solution constructed by Fife and Greenlee in [TJ5] (adapted to our choice of the function h) consists in 
adding an interior transition layer correction to expressions of the form ti + etj + e 2 tf, which approximate 
the solution u in the regions (notice that with our choice of the function h, we have Q± = f2+, f^2 = ^— , 
t\ = — 1 and ti = 1). This allowed Fife and Greenlee to construct an approximation U £ which yields an 
exact solution of (jTTJ) using a classical implicit function argument. No restriction on e are required, and 
the solution satisfies 

(8) u e — > —1 in f2_|_ and u e — > 1 in f2_ as e — > 0. 

Super-subsolutions were later used by Angenent, Mallet-Paret and Pclctier in the one dimensional case 
(see [7]) for construction and classification of stable solutions. Radial solutions were found variationally 
by Alikakos and Simpson in [5]. These results were extended by del Pino in [12j for general (even non 
smooth) interfaces in any dimension, and further constructions have been done recently by Dancer and 
Yan and Do Nascimento [16] . In particular, it was proved in [llj that solutions with the asymptotic 
behavior like ijHJ) are typically minimizer of J £ . Related results can be found in (TJ [2] . 

On the other hand, a solution exhibiting a transition layer in the opposite direction, namely 

(9) u e — * +1 in fi_|_ , m £ ^ -1 on fi_ as e — > 
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has been believed to exist for many years. Hale and Sakamoto [21] established the existence of this type 
of solution in the one-dimensional case, while this was done for the radial case in [13], see also [10] . The 
layer with the asymptotics in ([9]) in this scalar problem is meaningful in describing pattern-formation 
for reaction-diffusion systems such as Gierer-Meinhardt with saturation, see [T3] [18j [36] [41] [42] and the 
references therein. 

For one-dimensional or radial problems it is possible to use finite-dimensional reductions, which ba- 
sically consist in determining the location of the transition layer. In this kind of approach, the same 
technique works for both the asymptotic behaviors in |S]) and ©: the only difference is the sign of the 
small eigenvalue (of order e) arising from the approximate degeneracy of the equation (when we tilt the 
solutions perpendicularly to the interface). This makes the former solution stable and the latter unstable. 

On the other hand, one faces a dramatically different situation in higher-dimensional, non-symmetric 
cases. This is clearly seen already linearizing around a spherically symmetric solution of {T]) (with profile 
as in (O), as bifurcations of non-radial solutions along certain infinite discrete set of values for e — > 
take place, as established by Sakamoto in [42]. This reveals that the radial solution has Morse index 
which changes with e (precisely diverges as e — > 0, as shown in [U]). This poses a serious difficulty for 
a general construction. A phenomenon of this type was previously observed in the one-dimensional case 
by Alikakos, Bates and Fusco [3 in finding solutions with any prescribed Morse index. 

In [15] . del Pino, Kowalczyk and the third author considered the two-dimensional case, constructing 
transition layer solutions with asymptotics as in ([5]), while in this paper we extend that result to any 
dimension. Our main theorem is the following. 

Theorem 1.1 Let £1 be a smooth bounded domain o/lR™ (n > 2) and assume that a : O — * (— 1, 1) is a 
smooth function. Define K , f2 + and f2_ to be respectively the zero-set, the positive set and the negative 
set of a. Assume that Va 7^ on K and that K n dQ = 0. Then there exists a sequence Ej — ► such that 
problem ([3]) has a solution u ej which approach 1 in Q+ and — I in f2_. Precisely, parameterizing a point 
x near K by x = (jj, (), with y £ K and ( = d(x, K) (with sign, positive in f2+ ), u £j admits the following 
behavior _ 

u E] (y, ()=h(J- + <S>(y)^j + 0(e 3 ) as j -> +c^. 
Here $ is a smooth function defined on K and H (£) is the unique hetheroclinic solution of 
(10) H"+H-H 3 = 0, H(0)=0, H(±oa)=±l. 

As in [T3], [T5], [1B]-[1H]) [S] and other results for singularly perturbed (or geometric) problems, the 
existence is proved only along a sequence Ej — > (actually it can be obtained for e in a sequence of 
intervals (aj,bj) approaching zero, but not for any small e). This is caused by a resonance phenomenon 
we are going to discuss below, explaining the ideas of the proof. 

To describe the reasons which causes the main difficulty in proving Theorem ll.il we first scale problem 
^ using the change of variable x 1— > ex, so equation ([3]) becomes 

, . f Au=(n-(i(a))(ii 2 -l) inn e , 

1 ' I |; = on dn s , 

where f2 e = ~f2. Near the hypersurface K £ := \K, we can choose scaled coordinates (y, £) in Sl £ with 
y £ K e and £ = dist(x, K e ) (with sign), see Subsection 12.21 and we let u £ denote the scaling of u e to fl £ : 
with these notations we have that u e (y,() = u e (y,e() ~ H((). The function if(C) = H (dist(x, K e )) for 
x £ tt £ can be then considered as a first order approximate solution to pip , so it is natural to use local 
inversion arguments near this function in order to find true solutions. For this purpose it is necessary to 
understand the spectrum of the linearization of (jlip at approximate solutions. 

Letting L e be the linearization of (|11|) at u e , it turns out that L e admits a sequence of small positive 
eigenvalues of order e. Using asymptotic expansions (see Section [3] and in particular formula (|72|) ) . 

one can see that this family behaves qualitatively like e — e 2 Xj, where the Xj's are the eigenvalues of 

2 

the Laplace-Beltrami operator of K. By the Weyl's asymptotic formula, we have that Aj ~ j"- 1 as 
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j — » +00, therefore we have an increasing number of positive eigenvalues, many of which accumulate 
to zero and sometimes, depending on the value of e, we even have the presence of a kernel: this clearly 
causes difficulties if one wants to apply local inversion arguments. Notice that, by the above qualitative 
formula, the average spectral gap of resonant eigenvalues is of order e 1 ^ - . For the case n — 2 (considered 
in [IS]) this gap is relatively large, so it was possible to show invertibility using direct estimates on the 
eigenvalues. However in higher dimension this is not possible anymore, and one needs to apply different 
arguments. 

To overcome this problem, we use an approach introduced in [28], [29] (see also [25], [26], [27]) to 
handle similar resonance phenomena for another class of singularly perturbed equations. The main idea 
consists in looking at the eigenvalues (of the linearized problem) as functions of the parameter e, and 
estimate their derivatives with respect to e. This can be rigorously done employing a classical theorem 
due to T.Kato, see Proposition 13. 3[ and by characterizing the eigenfunctions corresponding to resonant 
modes. Using this result we get invertibility along a suitable sequence Ej — > 0, and the norm of the inverse 

n+l 

operator along this sequence has an upper bound of order 2 (consistently with the above heuristic 
evaluation of the spectral gaps). This loss of uniform bounds as j — » +00 should be expected, since more 
and more eigenvalues are accumulating near zero. However, we are able to deal with this further difficulty 
by choosing approximate solutions with a sufficiently high accuracy. 

Fixing an integer k > 1 and using the coordinates introduced after (111[) , from the fact that a vanishes 
on K, one can consider the Taylor expansion 



a(ey,eO=e(b(sy) + J2«) l h(£y)+Hy,0 with b(y,Q < C k \e(\ 



1=2 

and look at an approximate solution of the form 



ifc+i 



,{y, C) = H(( - $(ey)) + J2 ^{ey, C ~ Hey)), 



i=l 



for a smooth function &(ey) = <&o(ey) + Y^i=i £l ^i{ £ y) defined on K and some corrections hi defined on 
K x M+. Using similar Taylor expansions of the Laplacc-Beltrami operator in the above coordinates, see 
Subsection 12.21 the couple (hj, 3>j-i) for j > 1 can be determined via equations of the form 



(12) 



C h! = -K{ey)H'{s) + {s + $ )b(ey){l - H 2 (s)) 

C hj = Sj-xbiey)^ - H 2 (s)) + $ k (s, $ , • • • , ^-2,^1, ■ ■ ■ , fy-i, &i, ■ ■ ■ , bj), for j > 2, 



where Cqu = u" + (1 — 3H 2 )u, $k is a smooth function on its argument, and s — C — ®(£y)- (UU) is always 
solvable in hj by the Fredholm alternative if we choose properly the functions $1 . 

Such an accurate approximate solution allow us, using the above characterization of the spectrum 
of the linearized operator and the bound on its inverse, to apply the contraction mapping theorem to 
find true solutions. Specifically for the homogeneous Allen-Cahn equation, a related method was used 
in [3T] to study the effect of Oft on the structure of solutions to . Some common arguments are here 
simplified, and we believe our approach could also be used to handle general nonlincarities as in [19 1. 

The paper is organized in the following way: in Section [5] we collect some preliminary results concerning 
the profile H, we expand the Euclidean metric and the Laplace-Beltrami operator in suitable coordinates 
near K e , and recall some well-known spectral results. In Section [3] we first construct approximate so- 
lutions, and then derive some spectral properties of the linearized operator characterizing the resonant 
eigenfunctions: this is a crucial step to apply Kato's theorem. Finally, Section [4] is devoted to the proof 
of our main result. 
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2 Notation and preliminaries 



In this section we first collect some notation and conventions. Then, we list some properties of the 
hetheroclinic solution H, and we expand the metric and the Laplace-Beltrami operator in a local normal 
coordinates. Finally we recall some results in spectral theory like the Weyl asymptotic formula. 

Notation and convention 

We shall always use the convention that capital letters like A,B,... will vary between 1 and n, while 
indices like i, j, . . . will run between 1 and n— 1. We adopt the standard geometric convention of summing 
over repeated indices. 

j/i, ... , 3M-i will denote coordinates in R n , and they will also be written as y = (yi, . . . , y n -i), while 
coordinates in 1™ will be written x = (y, C) € M n_1 x R. 

The hypersurface K will be parameterized with local coordinates y = (y 1 , . . . ,y n _i). It will be convenient 
to define its dilation K e := -K which will be parameterized by coordinates (3/1, ... , y n -x) related to the 
3/' s simply by y = ey. 

Derivatives with respect to the variables y, y or £ will be denoted by dy, d y , and for brevity we shall 
sometimes use the notations di for d Vi . When dealing with functions depending just on the variable £ we 
will write H' , h' , ■ ■ ■ instead of OqH, OqIi, 

In a local system of coordinates, (<7jj)ij are the components of the metric on K naturally induced by 
R" _1 . Similarly, (jjAB)AB are the entries of the metric on fl in a neighborhood of the hypersurface K. 
( K i)ij wn l denote the components of the mean curvature operator of K into R n_1 . 

For a real positive variable r and an integer to, 0(r m ) (resp. o(r m )) will denote a function for which 
°^~ m - remains bounded (resp. °^ r m - tends to zero) when r tends to zero. For brevity, we might also 
write O(l) (resp. o(l)) for a quantity which stays bounded (resp. tends to zero) as e tends to zero. 

2.1 Some analytic properties of the hetheroclinic solution H 

In this subsection we collect some useful properties of the hetheroclinic solution H to (fTUjl . Note first 
that H can be explicitly determined by 

(13) H(C) = tanh(^C), 
and moreover the following estimates hold 

H(0 - 1 = -Aoe-^M + 0{e-^^) for C -» +00; 

(14) { H(() + l = A e-^+O(e-( 2 ^W) for C -> -00; 

H'(C) = V2A e-^\ + 0{e-^^) for |C| -» +00, 

where Aq > is a fixed constant. We have the following well-known result (we refer to Lemma 4.1 in [33] 
for the proof). 

Lemma 2.1 Consider the following eigenvalue problem 

(15) <p" + (1 - 3H 2 )<P = X(f>, <f> e If^R). 

Then, letting (Xj)j be the eigenvalues arranged in non increasing order (counted with multiplicity) and 
(4>j)j be the corresponding eigenf unctions, one has that 

(16) A x =0, <j)i^cH'; A 2 <0. 
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As a consequence (by Fredholm's alternative), given any function ip S L 2 (R) satisfying J R ipH' = 0, the 
following problem has a unique solution <fi 

(17) 4' + (1 - 3H 2 )cf> = ip in E, / H'(j> = 0. 

JR 

Furthermore, there exists a positive constant C such that ||0||//i(r) < CH^Hi^r). 
We collect next some useful formulas: first of all we notice that 

(18) H' = — (1 - H 2 ) and H" = -V2HH'. 

Moreover, setting 

(19) C u = u" + (1 - 3H 2 )u, 
we have that 

(20) Cq(HH') = -3V2H(H') 2 . 
2.2 Geometric background 

In this subsection we expand the coefficients of the metric in local normal coordinates. We then derive 
as a consequence an expansion for the Laplace-Beltrami operator. First of all, it is convenient to scale 
by - the coordinates in equation ([3]) to obtain 



(21) 



Au = (u — a(ex))(u 2 — 1) in E , 
|^ = on dn e 



where we have set f2 e — i£l. Following the same notation we also set K £ = ^K, and for 7 S (0, 1) we 
define 

S £ = {x e fl £ : dist(x,K e ) < e" 7 } ; I £ = [-e" 7 ,^ 7 ]. 

We parameterize elements x € S £ using their closest point y in K £ and their distance £ (with sign, positive 
in the dilation of f2 + ). Precisely, we choose a system of coordinates y on K, and denote by n(y) the 
(unique) unit normal vector to K (at the point with coordinates y) pointing towards f2_. Choosing also 
coordinates y on K £ such that y = ey, we define the diffeomorphism T £ : K £ x I £ — > S £ by 

(22) T £ (y,C) = y + Cn(ey). 

We let the upper-case indices A, B,C, . . . run from 1 to n, and the lower-case indices run from 1 

to Ti—l. Using some local coordinates (j/j)i=i,...,n— 1 on K e , and letting tp £ be the corresponding immersion 
into K™, we have 

f ff (f. = ^ + £ C^(ey) = ff(y) + eC^(ey)^f(y), for < = 1, . . . ,„ - 1; 
\^(y,C) = n( £2/ ). 

where are the coefficients of the mean-curvature operator on K. Let also (^)jj b e the coefficients 
of the metric on if e in the above coordinates y. Then, letting g — g £ denote the metric on Q £ induced 
by R™, we have 

(23) 9AB -^^-f^ ° 
where 

9ij 



dxA ' 9xb ) V 1 
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Note that also the inverse matrix {g AB } decomposes as 



g AB 



(g ij ) o 
o 1 



From the above decomposition of gAB (and g AB ) and for u defined on S e , one has 
A g u = g AB u A B + -j== d A (g Ai V det g) u b 

(24) = u a + g ij Uij + ^=5 C (v/d^) u c + ^=3 4 (g ij VfeUj) u 3 . 
We have, formally 

(25) detg = det(g- 1 g)detg= (detf) (l + eC tr (g^a)) + o(e), 
where 

= K t9kj + KjVii- 

There holds 

(g _1 a) is = 7/~ ''<>,< = V sj (^Vkj + K j9u) , 

and hence 

(26) tr (g-'a) = g ij {^9kj + = 2g ij n k g kj = 2 K *. 

We recall that the quantity n\ represents the mean curvature of K, and in particular it is independent of 
the choice of coordinates. 

We note that the metric gAB can be expressed in function of the metric g^, the operator Kj, and the 
variable e£. Hence, fixing an integer fc and using a Taylor expansion, we can write 

fc 

(27) -==d c ^r g = Y J ^C e - 1 Gt{ey) + G{ey,0, 

Vdet g ^ 

where Ge : K — ► K arc smooth functions, and G satisfies 

(28) G(-,C) ' < C ktm \C\ k e k+1 , (G/ e , 

where C k ^ m is a constant depending only on K, k, and m. Again (and in the following), when we write 
|| • ||' we keep the variable £ fixed. In particular, from the above computations it follows that 

(29) G^ey) = K(ey) := 4(ey). 

We need now a similar expansion for the operator A g : fixing the variable £ g I £ , the metric g{y,C) = 
g £ (y, C) induces a metric g £ ^ on K in the following way. Consider the homothety i e : K — ► K e . Wc define 

<?e,C = £ 2 *e.9e(->C), 

where i* denotes the pull-back operator. Basically, we are freezing the variable £ and letting y vary. 
Fixing an integer fc, for any smooth function v : K — > R we have the expansion below, which follows from 
(|24|) . reasoning as for (|27|) 

fc fc 

(30) A §e , { w = Y^eCfLtv + L E , k+1 v = A K v + J2«Y L ev + e k+1 L E , k+1 v. 

e=o i=i 
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Here L £ ^+i are linear second-order differential operators acting on y and satisfying 

(31) \\Uv\\' cm{K) < C m \\v\\' Cm+ . iK) ; L £ , k+1 v ' < C m \C\ k+1 \\v\\' Cm+2(K) 
for all smooth v, where C m is a constant depending only on K, k, and m. 

Consider now a function u : S e — > R of the form 

(32) u(y,C) = u(sy,C), V G K e , C S I e . 
Then, scaling in the variable y, we have 

A gMv, = u cc (ey, C) + —^=^ d < (V det 5) «f N, + ^^^(ey, 0- 
Using the expansions (f2"7| . (|3"0)) together with (j2"5)) , the last equation becomes 

A s ^(y, = u C( (sy, + (ek + J2 ^C^GiJ u c {ey, Q + G(ey, ()u c (ey, () 

k 

(33) + e 2 A K u + Y / £ 2+e C e L e u(ey,C)+£ k+3 L £ , k+1 u(eyX). 

t=i 

2.3 Spectral analysis 

We define the scaled Euler functional J e (u) in fl E by 

(34) J e {u) = \ ( |Vw| 2 + / F{ex,u)dx, with F(x,u) := [ (s - a{x))(s 2 - 1) ds. 
We set for brevity 

(35) b(y):=d n a(y,0) 

and we notice that by our choice of n, we have b > on K. Now, we let ipj and Xj be the eigenf unctions 
and the eigenvalues (with weight b) of 

-Ajcifj = Xjb(y)(fj. 

The Aj's can be obtained for example using the Rayleigh quotient: precisely if Mj denotes the family of 
j-dimensional subspaces of H 1 ^), then one has 



(36) Xj = mf sup f = sup mf f , 

where _L denotes the orthogonality with respect to the I? scalar product with weight b. We can estimate 
the Xj using a standard Weyl's asymptotic formula ([9]), one has 

(37) X 3 ~ C K ,bj^ as j -> +00, 
for some constant Ck,& depending only on K and 6. 

3 Asymptotic analysis 

This section is devoted to the construction of approximate solutions to (|2ip , and of approximate eigen- 
functions (and eigenvalues) in the £ component (see the coordinates introduced in (|22|) ) of the relative 
linearized equation. Then we characterize, via Fourier analysis, the profile of resonant eigenfunctions in 
both the variables y and (. 



3.1 Approximate solutions and eigenfunctions 

In this section, given any integer k > 1, we construct an approximate solution Uk, e to problem (|21|) . which 
solves the equation up to an error of order e k+1 . Using the above parametrization (y, £) in S e , we make 
the following ansatz 

k 

(38) « M (y,C) =H((-$(sy))+J2 £ih i( s y^-Hey)) m S e , 

»=i 

where is the hetheroclinic solution of (fTOf and where ^(ej/) = $o(ey) + e *^i( e J/) f° r some smooth 

functions $j defined on K. The corrections (/ij)j and are to be constructed recursively in the index 
i, depending on the Taylor expansion of a and the geometry of K. Since all the hi a will turn out to have 
an exponential decay in £, itfc i£ can be easily extended (via some cutoff functions) to an approximate 
solution in the whole f2 e , see (|102[) below. 

We first determine hi by solving the equation up to an error of order e 2 . To this aim, we expand the 
function a in powers of e as (notice that a(ey, 0) = 0) 

k 

(39) a(ey, e() = eb(ey)C + £>C)'Me») + Kv, *), 
where b{y 1 z) is smooth and satisfies 

\b{y,z)\<C k \ez\ k +\ 

Using the above expansion of the metric coefficients and the Laplace-Beltrami operator, see in particular 
(l33l) . setting s = C, — we obtain that the term (formally) of order e in the equation is identically zero 
if and only if the correction hi satisfies 

(40) C hi := (hi) ss + (1 - 3H{s) 2 )hi = -K(sy)H'(s) + (s + $ )b(ey)(l - H 2 (s)). 

By the asymptotics in (fT4|). the right-hand side is of class L 2 in R and, by Lemma \2. 11 (j40"|) is solvable 
provided the latter is orthogonal in L 2 to the function H'(s). Since H'(s) is even in s and since b(ey) > 0, 
this is possible choosing <j>o(ey) so that 

(41) b{ey)$a(ey) = K (ey) ,J~°° = ^K(ey). 

Moreover one can prove, using standard ODE estimates, that hi has the following (regularity properties 
and) decay at infinity 

(42) \d l s d™hi(£y,s)\<C m e m (l + \s\)e-^, 1 = 0,1,2, m= 0,1,2,..., 
where C m depends only on m, a and K. 

To obtain the other corrections ($i)i an d (hi)i one can proceed by induction, assuming that N > 2, that 
$o, . . . , $Ar-2 and hi, . . . , /i/v-i have been determined, and that (hi)i<N—i satisfy 

(43) \d l s d™hi{ay,s)\ <C m e m (l + \s\ di ) e -^, i < N - 1 2 = 0,1,2, m = 0,1,2,..., 

where C m depends only on m, a, K and di only on i. When we expand the equation (|21|) for u — un. £ 
in power series of e, the couple (/ijv, ^jv-i) can be found reasoning as for (/ii,$o) : indeed, considering 
the coefficient of e N in this expansion, one can easily see that /ijv satisfies an equation on the form 

(44) £ h N = $ A r_ 1 fe(ey)(l - H 2 (s)) + $ N (s, $ , ■ ■ ■ , $jv-2, • • • , h N - X ,b x , b N ), 
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where 3jv is a smooth function on its argument. Reasoning as for hi this equation is solvable provided 
the right-hand side is L 2 -orthogonal to the function H'(s). This is indeed true choosing $at_i so that 

- f + °° H'(s)3 N (8,-)d S 
b(Ey)9 N -i(ey) = —r^ . 

Furthermore, one can show that /ijv satisfies regularity and decay estimates as in (|43|) . Reasoning as 
in Section 3 of [29j one can check that the above formal estimates can be made rigorous, and that the 
exponential decay of the corrections yields the following result. 

Proposition 3.1 Given any integer k > 1 there exist a function u k ^ £ : S £ — > M which solves equation 
(I2ip up to an error of order e k+1 . Precisely, setting 



(45) & e (u) = Au- {u-a(ex))(u 2 - 1), 
there exist a polynomial P k (C) such that 

(46) \6 e (u kie (ey,C))\<e k+1 P k (C)e-^ lcl in S e . 
Moreover, the following estimate holds 

(47) \d l s d™u k , e (s y> 0\ <C m e m P k (C)e-^, 1 = 0,1,2, m = 0,1,2,..., 
where C m is a constant depending only on m, a and K . 

We will look at solutions u of (|2~Tj) as small corrections of u k:£ (suitably extended to fl e via some cutoffs 
in £, see (I102p below), namely of the form 

u = u kyE + w, 

for w small in a sense to be specified later. The equation 6 £ (u) = is then equivalent to 

L e (w) + e e (u k , e )+K(w)=0, 
where L e is nothing but the linearized operator at the approximate solution u k ^ e 

(48) L E w := Ag F w + (1 — 3ul £ )w + 2a(ex)u k:£ w, 
and where Af e is the remainder given by nonlinear terms in 6 £ , namely 

(49) Af E (w) := -(3it fe , E - a(ex))w 2 - w 3 . 
It is also convenient to define the following linear operator 

(50) L e w := w cc + J_ d c (^det g] w,; + (1 - 3m 2 £ )w + 2a(ey, ()u ky£ w. 

Vdetg v / 

In particular, using the expansion (|33p . the operators L £ and L £ are related by the following formula: 
setting w(y, £) = w(ey, C,) one has 

(51) L £ w = h £ w + e 2 A K w + e 3 £ 3 , e w, 

where £3^ consists of the last two terms in (|33l) (replacing u with w). Precisely, £3^ is a linear differential 
operator of second order acting on the variables y, which for every integer m satisfies 

(52) £ 3 , e v < P(0\H 

G (K) 
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Here P(Q is a polynomial in £ with fixed degree, and coefficients depending only on m. 

We want next to derive some formal estimates on the following eigenvalue problem 

(53) h £ v = fxv in I e , 

with zero Dirichlet boundary conditions. It follows from Lemma 12.11 that the eigenvalues either stay 
bounded away from zero, or converge to zero as s — > 0: we are interested in the latter case. We argue 
heuristically expanding (|53[) at first order in e. In the limit e — > 0, we have /i = with corresponding 
eigenfunction H' , therefore it is natural to look for approximate eigenfunctions of the form 

(54) V = H' + eH 1 , 

and eigenvalues /x = e~p + o(e). We impose that Hi is orthogonal to H' in L 2 (K). Therefore the 
approximate eigenvalue equation formally becomes 

(55) C H l = -2b(ey){s + <5> )HH' - H"n(ey) + QHH'hi + JIH' + o(l). 

As for (|4"0)l . solvability is guaranteed provided the right-hand side is orthogonal in L 2 to H' . Using 
the oddness of H, formulas (|18[) . ([2(7]) . (|4"0)) and the self-adjointness of Cq we find that orthogonality is 
equivalent to 

, . LsH(s)(H'(s)) 2 ds n , 

With this choice of /z, the function Hi is defined as the unique solution of 

C Hi = -2b(ey){s + $ a )HH' - H" n{ey) + QHH'hi + V2b(ey)H'. 

From the exponential decay of H' and hi (see (j42|) ) we deduce that Hi satisfies estimates similar to (j43|) . 
Using this fact and the rigorous expansions in (|27)), (|3^|) . we then derive the following estimate 

(56) !, e (^)=eJlH' + e 2 R E (ey,C), 
where the error term R e satisfies 

(57) |i? £ (ey,C)| <P(C)e-^ ICI 

for some polynomial P(Q- Also, from the regularity and the decay of Hi we have 

(58) \d l s d™^{ey,s)\<P{QE m+1 e-^, 1 = 0,1,2, m = 0,1,2,. ... 
3.2 Characterization of resonant eigenfunctions 



We characterize next the eigenfunctions of L e , see (|48|) . corresponding to small eigenvalues. Let us recall 
first the definition of (<fj)j and (Xj)j in Subsection 12.31 

Lemma 3.2 Let X e = 0(e 2 ) be an eigenvalue of the linearized operator L £ in S £ with eigenfunction <f> 
and weight b, namely 

L £ <j> = X £ b<p in S £ 

(and with zero Dirichlet boundary conditions). Let us write the eigenfunction </> as 

<f>(y,Q = <p(ey)*(£y,0 + <t> ± {y,0, 

with ^ defined in (|54|) and with (f> satisfying the following orthogonality condition ( we are freezing the 
y variables in the volume element) 



(59) / *(ey,C)<^(y,CWe(C) = for every y e K £ 

Then one has \\4> ± \\h 1 (s e ) = o(e) |] 0|| j^i ^ as e tends to zero. 
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Proof. We notice first that, since \& = H ' + o(l) in H 1 (R), the operator Co is negative definite on (p 1 - 
by Lemma 12.11 Therefore, using the estimates on the metric g £ in Section [21 we find easily that there 
exist a constant C > such that 



(60) 



1 



^(y,OL £ <t> ± (y,OdV gB (y,0 < -^ll^ll^s.) 



Let us write the eigenvalue equation L £ <j) = \ £ bcj> as 

i e 0-L = -L e {<p%) + KH 1 - + KbtpV. 
Multiplying by (p- 1 , integrating over S e and using (I59|) we obtain 

/ <^(y,C)£ e <^(y,C)^ = - / ^(j/,C)i e ^N)*N,C))^+A E / fe(£y)^(2/,C) 2 ^ 

By ([60)1 (and the smallness of A e ) it then follows 



(61) 



Mlffi(S e ) - C 



^{y,£)L e {v{ey)V{ey,Q)dV gi 



Now by (|51[) and (f56"| we can write L e ((/? X I') as 

L £ (^*) = L e (^)+e 2 A K (^)+e 3 £ 3!£ (^) 
(62) = <p (eV2b(ey)H' + e 2 R £ (ey, C)) + e 2 A K (^) + £ 3 £ 3 ,e(^)- 

Then, again by (|59p . we have 



± ( 2 /,C)ie(^N)*N,C)) 



(63) 



+ 



^(y,0&K&(ey)*(Ey,0) 



(e 2 R £ (sy,CM£y) + £ 3 hA<p(zy)y(zyX))) ^ 



where R e is as in (|5T|) . We first estimate the second term: since 2,%^ is a second order operator in y 
satisfying the bound (|52p , by an integration by parts we find 



<Ce< I P{Q\V v {ip{ay)*[eyX))\\Vy^\dV a M£)- 

K E xI s 



Therefore, using the Holder inequality, the change of variables y = ey, (|58|) and the estimate on R e we 
find 



(64) 



e 2 R e (ey,(Mey) + e 3 £ 3 , e (v(ey)*(ey,())) ^ dV gc (y,C) 



£ 2 



for some positive constant C. It remains to estimate the first term in (|63|) . To this aim we decompose (p 
as (see the above notation) 



(65) 



£2/) = J2a j tp J (ey), 



for some real numbers a,-. One can write 



<^(y, 0A K (<p(ey)9(ey, C)) dV ge (y, C) 



<^(y, C)*(ey, C)Ajf ^(ey) dV 9e (y, C) 



+ / <j> ± (y,CM£y)A K ^(ey,C)dV gs (yX) 
Js E 

+ 2 [ ^{y,QV K ^{ey,C)-V K ip{ey)dV gB {y,0- 
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The first term vanishes by (|59j) . Hence, using the Holder inequality, (|58|) and reasoning as for (|64|) we 
obtain 

1 



<j> ± (yX)A K (<p(ey)y(y,())dV a Ay,0 
which by estimates (|63| and (|64|) implies 
(66) 



< C—^rW^WHHs^WvWmiK), 

£ 2 



^(y,C)i £ (^(y)*(y,C))^ e (2/,C) 



< C —^rW^WH^s^WWmiK)- 

e 2 



Using then (fBTj) and the latter equation together with the Weyl's asymptotic formula (see Subsection 
2.3|) , we then obtain 

C 



(67) 



Now, we rewrite the eigenvalue equation as 

(eV2bip + e 2 A K ip^ * = + A £ V* - s 2 ipA K ^ - 2e 2 V K ¥> ■ Vk* 

We use again the above decomposition <p{ey) = ( £z ); we define the integer j £ (depending on 

e) to be the first j such that e 2 Xj > £=. We multiply this time the last equation by J2j>j a 3 l Pj^ 
and we integrate over S £ . By the self-adjointness of L £ , (|58[) . and a similar argument as for 
incorporating the term involving £3.^ into the left-hand side we obtain that 

(i + o(*))i E ^ * c ^ + w) {^=1 E a u \J=t E < 



(68) 



The last term can be estimated reasoning as for (j6"7| : since Aj ^> 1 for j > j e , we find that 



JS * i>i. 



< Ce\\^\\ HHSs) UfJ^j 



From the above estimates and the fact that A E = 0(e 2 ) we obtain 



(69) 

Then, writing ipty as 



j>je 



* E • E "..^v 1 '- 



using (|6T)l and (155)) we find 

H^IUiCS.) < C(e 2 +ei)\\^\\LHs c )+Ce [ -L ^ £ 2 A, a j 



< CeiH^IU-cs.) + C £ 2 (||^|| i2(Se) + ||^|| ffl(Se) ). 
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This implies H^H^i^) < Cel\\(p9\\ L 2^ Se) , and noticing that ||^|||2 (Se) = H^ll!^) + \\y \\ L 2( Se 
achieve the desired estimate. ■ 

Our next task is to estimate the derivatives of small eigenvalues of the linearized operator L £ with respect 
to the parameter e. This will allow us to obtain invertibility of L e for a suitable family of small e. The 
prove of the main result can be then obtained by a direct application of the contraction mapping theorem. 
Using a result by T.Kato, see [25], page 445, which can be applied by the symmetry of L £ and elliptic 
regularity results (these ensure that the eigenvalues of L e are stable and semi-simple, according to the 
definitions in 22J), we have the following proposition. 

Proposition 3.3 For k £ N, k > 1, let Uk, £ be given by Proposition 1 3. 1\ and let L £ be defined in (|48|) . 
Then the eigenvalues X e of the problem 



(70) 



L £ u — X e bu in S e ; 
u = on dS E , 



are differ entiable with respect to e, and they satisfy the following estimates 
(71) Tl, E < ^ < Tl e . 

Is gM 2 ~ § u k,eVk,eU 2 ) dV gt 



Here we have set 



TL = inf 



ueH x ,u=£o J s bu 2 dVg e 

Is c ("I l V 9. M | 2 ~ &U k , e V Ke U 2 ) dV gt 



ueH x ,u^o J s ^bu 2 dVg E 

Vk, e = ; while H\ stands for the eigenspace for (|70p corresponding to the eigenvalue X. 

We next give a further characterization of some eigenfunctions of L e , in addition to the ones in Lemma 
concerning in particular the function ip. 



Lemma 3.4 Suppose the assumptions of Lemma \3.S\ hold true. Then, normalizing (f> by \\4>\\h 1 (s e ) = 1; 
decomposing ip as in (|65p and setting 

(72) Aj, e = V2e-e 2 Xj, 
as e — > we have that 

\\ 3 ,e\>ei |Aj, s |>e3 

Proof. We define the sets 

M,e = {j e N : Xj, e > e^} ; A 2<e = {j£N : Xj, £ <-£*}, 

and the functions 

<t>i = ^i(ey)*(ey, 0; <p2 = <p 2 (ey)V(ey, ()■ 

As one can easily see, from the estimates on g E in Subsection 12.21 and from the decay of see ([S"5]) . as 
e — > there holds 

(73) \m& H s e ) = ^r-(ca j K V 2 + C x /^IV^I 2 ), 
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where Co = J M (H') 2 + (H") 2 and C\ = J R (H') 2 . Similar formulas hold true for <f>i and (j>2, and hence 
these two functions stay uniformly bounded in H 1 (S e ) as e tends to zero. 

We multiply next the equation in (|70j) by 4n, I = 1,2: from the orthogonality of lf 1 and Tp 2 on K (with 
weight 6), (|59|). an integration by parts and the above arguments we get 



0(e 2 )\\Ul^) + 0(s 3 )\\vnLHSc)Ui\\LHs c 
Using (fM]) (replacing ip by Tp{) we deduce 



iL e (f>dV gi 



(^ + 0- L )i e ^. 



< C- 



£ 2 



for Z = 1.2. 



Also, from the expression of L £ (see (|62]l ). the subsequent estimates and some straightforward computa- 
tions one finds 



e<Pl 



1 + 0(1) 



m E * e 



Then from the last three formulas, Lemma 13.21 and the normalization on <b we obtain 



(74) 



lj E «^ = ^ 2 ) [-L E «i] +0(e) * E ^ 



for Z = 1,2. 

We next further split the sets A^ e as Ai :E = Ai >£ U Ai :E , where 

Ai. e = {jeA le : |£ 2 Aj|<e*}; A ltS = Ai >e \A hs , 



so from ([74")) we have clearly 

-sn E a i A ^ 



(75) 



O(e) f ^ E 



Obviously, since Xj t£ — \/2e — £ 2 Aj, for j S -A; j£ the ratio ^— ^ stays uniformly bounded from above and 

below by positive constants. Therefore, using the elementary inequality \xy\ < 8\x\ 2 + j\y\ 2 with 8 small 
and fixed, we can absorb the last term into the left-hand side of the latter formula, obtaining an error of 
the form 0(e 2 /6). Therefore, using also the definition of we find 

* E d$\te=0(e*)(l E o 2 ] + 0(e 2 ). 



jeAi 



jeAi 
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By our normalization on <j>, see also the comments after (|T3[) . the argument inside the last bracket is 
uniformly bounded as s — » 0, and hence we obtain the second assertion of the lemma. 

To obtain also the first one we notice that |Aj, e | > ei for j S A^ e , so we find 

as e — » 0. This concludes the proof. ■ 

Now, using the above lemma, we can estimate the derivatives of small eigenvalues of L e with respect 
to e. Precisely, we have the following result. 



Lemma 3.5 Let X be as in Lemma \3.SX Then, for e sufficiently small A is differentiable with respect to 
e, and there exists a negative constant CK,b, depending only on K and b, such that its derivative (which 
is possibly a multi-valued function) satisfies 



d\ 

o c K,b 

OS 



o(l) as e — > 0. 



Proof. The proof is based on Lemma [3.21 Proposition 13.31 and Lemma [3.41 Since we want to apply 
formula (|71[) (in our previous notation) to the function 

w = 0= a iVA £z ) \V + <l>- L =<p(H' + eH 1 ) + <l)- L , 

\\^-,s\>ei ) 

we need to estimate the two quantities 

(76) J ^ (^-~ £ \y gM 2 ' ^k,eV k ,eU 2 ^j dV ge - J^bu 2 dV ge . 

Here the function Vk, e is defined as Vk, e = dU g £ E (£■), where uu,e '■ zS e — *■ R is the scaling Uk. e {x) — Uk. £ (£x). 
We claim that, normalizing u with ||ti||fl-i(s s ) = 1 (this condition was required in Lemma [3.4|) . the following 
estimates hold 

(77) j s (-- £ \V g u\ 2 ~ 6u k , £ v K£ u 2 ^j dV ge = ^ b(y)<p 2 + o(l); 



(78) / bu 2 dV ge = -^ [ b(y)^ + o(l) 

as e — ► 0, where c < and where C\ is defined after (f73|) . This together with (|7ip would conclude the 
proof of the lemma. 

Proof of ([77)1 and ([75)1 . First of all, recall that by our normalization and by Lemma 13.21 we have that 
ll < / ) ~ L ||-ff 1 (s e ) — °{ £ )- Therefore, using the expansions for g e in Subsection 12. 2[ some integration by parts, 
([58)) and the estimates in Subsection 13.11 one finds 

(79) - / ip{H' + effi) [p{H'" + eH'l) + en V {H" + eH[) + e 2 A K (p(H' + eHi)] (1 + e<) 

(H' 2 + 2eH'H 1 )(1 + en() + o(l). 



ip 2 (H + sh 1 ) 
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Since the arguments of H,H',hi and H\ are all translated by $o m Cj with the change of variables 
s = C — $0 and some elementary estimates (which use the exponential decay of H' , hi and Hi) we find 

/ (--\VgA 2 - 6u k , £ v k . £ uA dV gs = - f (p 2 (2H'H"' + 6sHH' 3 ) 

JS e V £ J e JJfjXl 

(80) +4/ ip 2 Hi (H'" + 3sHH' 2 ) + 6 f tp 2 ( shi(H') 3 + sI^HiH') 2 ) 

JK c xR JK e xR V / 

+ 2 f eipA K ipH l2 + 2 [ k(s + ® Q )H'H'"ip 2 + 6 f s(s + $ )KHH' a cp 2 + o(l). 

JK e xR JK e xR JK E xR 

In the latter formula all the arguments now are simply in s, with no more translation. Using equation 
(|10p. the oddness of together with some integration by parts, it is easy to see that the term of order 
- in the above expression is identically equal to zero. Let us consider now the terms of order which 
involve H\. Using the self-adjointness of the operator Cq and the following elementary identity 

£0 f — \=sHH') = H'" + 3sHH' 2 , 



V2 



we can write that 



if <p 2 Hi (H'" + ZsHH' 2 ) = — — f Lp 2 sHH'CoHi. 

J K e x~R v2 Jjf e xl 

Similarly, the terms of order which involve hi can be written (up to some integration by parts in the 
variable s) as 

6 f ip 2 ( shi(H') 3 + sh'iH(H') 2 ) = -6 / tp 2 ( 2sHH'H" + H (H') 2 )hi 
(81) = -6 \ (p 2 H{H') 2 hi + 12V2 [ ip 2 s(HH r ) 2 h 1 



2 
71 



I K E xR JK E 

[ ip 2 HH'£ (hi) + ^ [ ip 2 s(HH') 2 hi. 

JK s xR \2 Jk e xR 



Here again we have used the self-adjointness of the operator Co and the identities 

H" = -V2HH' and C ( ^=HH'\ = -6HH' 2 



.V2 

Regrouping the above terms, and using the oddness of H one finds 
J ^~^|V 9e w| 2 - Quk^Vk.eU 2 ^ dV 9s 
= ~^=J V 2 sHH' (c Q Hi - QHH'h^j 
+ 4=/ <p 2 HH'£ (hi) + 2 f e<pA K tpH' 2 

(82) + 2 / k(s + <S> Q )H'H'"(p 2 + 6 / s{s + <S> Q )KHH' 3 (p 2 + o(l) 

J K E xR J K € xR 

[ if 2 sHH'(c Hi-6HH'hi) + [ ip 2 HH'C {hi) 

JK c xR V / \/2 JK e xR 



4 

71. 



+ 2 [ eipA K ipH l2 + 2 [ K<P ip 2 (3sHH 13 + H'H 1 ") + o(l). 

JK F xR JKrXM 
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Now, using equations (|40|) and (|55| we arrive at 

J (^"7l V ^ M l 2 _ 6u k , £ v k ^u 2 ^J dV gi 



/ <f 2 sHH , [2b(ey)(s + <f> )HH' + H"n(ey) - V2b(ey)H / 

J if e XR V 



V2 

'K e xl 



V2 



K (ey)H' + (s + <S> )b(ey)(l - H 2 )^J 



+ 2 eipA K ipH' 2 + 2 ^^(ZsHH 13 + H'H"')+o{l). 

JK s xM JK c xR 

Using again the fact that H is odd, by the vanishing of the last integral (as one can easily check) we 
obtain 

j s ^-^ g u\ 2 ~^M,eU 2 ^dV ge = -^jf if 2 sHH'^2b{ey)sHH' - V2b(ey)H'^j 

+ 4=/ ip 2 HH'sb{ey)(l- H 2 ) + 2 [ EtpA K ipH' 2 

+ 0(1). 

By an explicit computation of the integral we find 
so by Lemma [33] and some easy estimates, we find 



bLp 2 +o{\). 



Then we obtain (JS) taking c = (^tt 2 - ^) < 0. To prove (|78|) it is sufficient to use Lemma 13.21 the 
estimates on g £ in Subsection 12.21 and the decay of Hi , see (f54| . ■ 



4 Proof of Theorem 11.11 

In this section we first prove the invertibility of the linearized operator L £ using Lemma [33] and choosing 
carefully the parameter e. Below, Hq(S £ ) stands for the functions in H 2 (S £ ) with null trace on dS £ . 

Proposition 4.1 Let k > 1, let u k ^ £ be the approximate solution defined in Provosition \S.l[ and let L £ 
be the linearized operator at u k}£ , see (|48|) . Then there exist a sequence Ej — » such that L £j : H 2 (S £j ) — > 
L 2 (S £j ) is invertible and its inverse L~ : L 2 (S £j ) — > HQ(S Sj ) satisfies 



L- 1 



< Cs , for all j e N. 

£(^(S e );i?g(S e )) 3 



Proof. The proof is similar in spirit to the one of Proposition 4.5 in [5S]. As we will see, in order 
to study the spectral gap of L £ it suffices to find an asymptotic estimate on the number N £ of positive 
eigenvalues of L £ and to apply then Lemma 13.51 We denote by Ai, £ > A2, e > • • • > Aj, e > • • • the 
eigenvalues of L e , counted with multiplicity. The j-th eigenvalue Xj^ £ can be estimated using the classical 
Courant-Fisher formula 

J s uL £ udV g ~ f s uL £ udVg e 
(83) X jt£ = sup inf / a ; X jt e = w mf sup — 



m g Mj u^M.u^o J bu 2 dVg F ' 3 ' e MeMj-i u ±m,u^o J s bu 2 dV 3l 
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Here Mi represents the family of /-dimensional subspaces of Hq(S £ ), and the symbol _L denotes orthog- 
onality with respect to the L 2 scalar product with weight b. Notice that the inf and sup are reversed 
compared to (|36|) since the principal part of the operator has the opposite sign. 

We can find a lower bound of N e using the first formula in (|83| . Indeed, given a fixed S > 0, let j e be 
the largest integer j for which Aj j£ > Se. From ([57)1 and ([72)) we find that 

n — l 

(84) Aa (^r)' as " ' 

We can take a test function </> like ip^ with ^ = YjlLi a iVi- Actually, since we want to work in the space 
Hq(S e ) we need to add a suitable cutoff function in £. However, by the exponential decay of 'J, see ([58)) . 
these will generate error terms exponentially small in e. Therefore, for convenience of the exposition, we 
will omit these corrections. 
By (162)) we have that 

L e (/> = (eV2bip + s 2 A K <p)y + e 2 LpA K V + 2e 2 {V K (p, V K ^) + e 2 R e tp + £ 3 £ 3 , e (¥>#), 

where £3.^ and R e satisfy respectively the estimates ([5"2")) and (1571) . Reasoning as for ([51)) . ([55)) we find 

(85) / ^L^dV^^-^^iil + oim^ + Ois^ah I b<?dV ge = i±£^£ a a. 

JSe £ I JS C £ 

Defining M = span{tfi^> , I < j £ }, by the first formula in ([55)) and our choice of j e we have that 

Is. uL eUdV gF 



X Je . e > inf ? 7 " > 0, ase^O. 
^£M,u^Q J s bu 2 dVg e 



From (|84p and the last formula we then find the following lower bound 

n-l 

(86) N s >(l + o(l))l^—^\ zse^Q. 

\ CK,b£ J 

A similar upper bound can be obtained using the second formula in ([83)) : again, given a fixed S > 0, let 
j e be the smallest integer j for which Aj ;E < —Se. Still from ([37)) and (17211 it follows that 



Now let </> G Hq(S € ) be an arbitrary function orthogonal to M := span{Lpi^> , I < j e — 1}, and let us 



write it in the form cf> = ipty + <p with (f> as in Lemma 13.21 

We write as before tp = cti<fi, and split it as sum of the following two functions 



Using the second formula in (|83p we have that 

Jo uL e udV ge 



< sup , , . 

I »> e OU £ dV„ 

By the definition of (j) 1 - , the expansions of the metric g e in Subsection 12.21 and ([58)) one finds 

\\<t>\\h{s.) = y*\\h(s s ) + ll^lll=(s e ) = C 1 + o(l))(||^ 1 *||i2 (Se) + ||? 2 *lll»(s.)) + ll^lll^s.) 

2 

L 2 {S E 



(89) = (1 + 0(1))-^ J2 «I+(l + o(l))^rE^ + ll^" ll? - a '- ' as . - i!. 
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Using also the estimates in the proof of Lemma 13.41 multiplying <p by 6(e?/)</7 1 v E', using the orthogonality 
to M and integrating one can easily prove that 

(90) \WML*(s e )<o{l)\\<j>\\ L *{Se), ase^O. 

This, together with (|^|) . the fact that |Aj- e | < Ce for j < j £ — 1 and some easy computations imply 

(91) f {7p^)L s (Tp^)dV ge = o(e)\\ct>\\l 2{Se) . 
Now, by the self-adjointness of L e we can write 

j 4>L £ cbdV gs = j {Tp^>)L £ {Tp^>)dV ge + I (^)L £ {lp 2 ^)dV ge + [ ^L £ ^dV 3c 

(92) + 2 I {^)L e {rp 2 ^)dV gc + 2 f (<p 1 *)L e <t> ± dV ge +2 f (Tp^L^dV^. 
Let us first estimate the last two terms: from we obtain 



\!>3* 



(93) (lp 2 *)L e <p x dV ge < — W^W^sJ^W^K) < __||0-L|| H1(Se) ^(i + Ai)c 

JS S £ 2 £2 \ . 

Similarly, using the fact that |Aj| < y for I < j e — 1 we have 

II^iIIhi(a') < E Xia ^\ < — II^iIU 2 (-k:), 



\l<je~l 



and hence 



f Ce 2 Cei 

(94) / {^)L £ ^dV gc < -^tW^WhusjWiWruk) < -^W^WhhsJViLhk)- 

JS C £ 2 £ 2 

On the other hand, a similar argument as for the first formula in (|85[) yields 



(95) / @ 3 *)L s faiti)dV ge + o(l))X LE + 0(e 2 )] 

Moreover by the negative definiteness of L £ on cf) 1 - , sec (l60l) , we have that 

A-L r jJ-Jt^ ^- /"(— 1 II A-L II 2 



a?. 



(96) / ^L £ ^dV ge <-C-^\\^ Se) 



for some fixed constant C. It remains to estimate the term J s ((p 1 *S>)L £ (<p 2 ' ( &)dV ge . Using again (jl?2")) . 
(|58[) . the fact that J„ bifiipj = for i < j s — 1, j > j e , and some integration by parts we get 



(97) ^r- / fa¥)L B C¥> 2 tf)(iVk < £ 3 / |V^| |V^ 2 | + e 2 f (|V^| fo| + |^| |V? 2 | + |^| |%|) 

< £ 3 ( e W] fcJ +e 2 ( E « 2 ) 

\i<J S -l / V>Js J \l<3e-l J \l>3e 

+ - 2 (e (e«N +^ 2 (e « 2 ) (e« 2 
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We claim next that the terms in the right-hand side of (|92|) can be combined yielding 



C U^E A 



(98) J s <j>L s ct>dV gs — | — - \.:<<l !l--!|'7/.,s. , J -"(-.O'ill-^., 

To prove this, we show that the main terms in (f9"2"| are the ones given in (|95p . while all the 

others, listed in the left-hand sides of formulas ([50]) , (f9"3"| . (IM|) . ([57)) can be absorbed into the formers by 
the elementary inequality \ab\ < a 2 + b 2 . For example, for any small constant (3 (independent of e) we 
can write 

^r\\^\\m(s,) ( + A ')«M ^ C^ll W> + ~ ^ E ^ + A <W- 

Taking /? sufficiently small, and noticing that e 4 A; > e 2 Xi. £ + 0(e 3 ), from (|89[) we deduce our claim for 
this term. For the others, one reasons similarly, taking also ([90]) and the choice of j e into account. Now 
([89]) . ([90]) . ([98]) and again the choice of j E imply 

J J>L e 4,dV ge < --^ E «? - ^II^IIW) < 0. 
Therefore, by (|5Tj) we find the following upper bound on N e 

71-1 

iV £ <(l + (l))(^tij ase^O. 



Since 8 is arbitrary, the last estimate and (|86[) imply 

(99) 7V £ ~ Ci^e"^ as e -► 0, 

where we have set C^if = (<9~~) ' ■ 

Next, for / g N, let e ; = 2"'. From ((Ml) we get 

(100) N El+1 - N ei ~ Ci lJC - 2^) = Ci,k(2^ _ 1) £ "^. 

On the other hand it follows from Lemma 13.51 that the eigenvalues of L e which are bounded (in absolute 
value) by 0(e 2 ) are decreasing in e. Equivalently, by the last equation, the number of eigenvalues which 

_ n — l 

become positive, when e decreases from si to £;+i, is of order e ; 2 . Now we define 
Ai = {e e (ei +1 ,ei) : fcerZ £ ^0}; Bi = (e/+i, £;) \ Ai. 

By (llOOp and the monotonicity (in e) of the small eigenvalues, we deduce that card(.4;) < Ce l 2 , and 
hence there exists an interval (a;, b{) such that 



(ioi) K6 z )cB i; \h-at\yc- 1 ,r;t >c-^ 



_ 1 meas(S;) . - — l 

card(^4;) 



From Lemma 13.51 we deduce that L a l +b l is invertible and 

2 

c 



r-l 



< 
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This concludes the proof taking e, = 3 . 



Proposition [47T] gives us a localized version of the invertibility result we need. To have a global one in the 
whole domain £l £l we define a smooth cutoff function \e by 

' X e(t) = 1, for t < 
XS) = 0, for * > fe-T 
|x' e | < Ce^ and < Ce 2 ^. 



We next set 
(102) 



i + x £ (0K. £ (y,C)"i), infi+ 
-l - xe(0 (ukAv, + 1) , in K n \ n+. 



Now by (fl"4"]) and (|47[) we know that |tifc e | is exponentially close to 1 and its derivative are exponen- 
tially small for \(\ sufficiently large. This and (|4"6"]) imply that \\& £ ({ik, £ )\\L 2 (n e ) < Ce k+1 ~ and 
||S £ (ufe !e )|| L oo( 0e ) < Ce k+1 , see (|4"5"1) . where C depends only on K, b and k. 

Let us denote by L £ the linearized operator at Uk,e m ^e- Given a smooth positive extension 6 of b 
to r2, we consider next the eigenvalue problem 



(103) 



L e u — Xb(e-)u in fl £ ; 
#i = on dfl £ 



and we denote its eigenvalues by (Aj, E )j, counted in decreasing order with their multiplicity. 
By (|48|) . asymptotically away from K £ , the eigenfunctions u satisfy 

Am - (2 ± 2a(ea;) - \b)u = in (fi±) e , 

where (f2±) e stands for the ^ dilation of f2±. Since we are assuming a G (— 1, 1) in O, if A is bounded 
from below by — e, the coefficient of u in the above equation is negative. Hence, reasoning as in [28] , 
Lemma 5.1, one can prove that u has an exponential decay away from K e . 

Moreover, an argument based on the Courant-Fisher method, see Proposition 5.6 in [28], shows that 
there exists a constant C depending only on fl, K, a and 7 such that 



provided A,-. £ > —e or A,-. £ > — e. 



\X jte -X jtS \<Ce-r 

Here Aj )E are the eigenvalues of L e in S e , see the proof of Proposition [O] 

This and Proposition 14.11 allow us to prove the following result which guarantee the invertibility of 
the linearized operator for the range of the parameter e constructed above. 

Corollary 4.1 Fix k € N and let Uk,e> L £ be as above. Define H^(fl E ) to be the subset of H 2 (Q e ) consist- 
ing of the functions with zero normal derivative at dVL £ . Then for a suitable sequence Sj — » 0, the operator 

L Sj : H^(Q ej ) — > L 2 (fl £j ) is invertible and the inverse operator satisfies 
for all j e N. 



L: 



< 



£(L=(fi..);ffJ(fi,.)) 



Using the above results, we are in position to prove our main result, Theorem ll.il 

Proof of Theorem 11.11 Let (ej)j be as in Corollary 14.11 We look for a solution u e of the equation 
®e( u ) = of the form 

u e =u k>e + w, wGH 2 (n £ ). 
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For e = e h define the function F e : H*(Cl e ) H i°°(^ £ ) -> i?2( fi e) n L 00 ^) by 

(104) F £ H := -L- 1 [6 B (fi*, e ) - (3u fc , £ - a) u; 2 - u; 3 ] . 
We have that 

(105) 6 e (ufc,e + w) = F e (u;)=u>. 

We want to prove that F e is a contraction in some closed ball of iJ^(f2 e ) n L°°(f2 E ). We first define the 
norm ||| • ||| as |||u>||| = Hi^ll H"2(n e ) + ||w||l»(S) e )- Then, for r > 0, we introduce the set 

B r = {weH*(n e )nL°°(n s ) : \\\w\\\<r}. 

Applying a standard elliptic regularity theorem and using Corollary 14.11 one can prove that there exists 
positive constants C (depending on fi, K and a) and d (depending on the dimension n) such that 

(106) III^MIII < Ce~ d (e k+1 -^ + ||M|| 2 ) ; 
\\\F e { Wl )-F e {w 2 )\\\ < Ce-'dlKIH + IIIHIDdlK-HII). 

for e = Sj and w, wi,W2 € -ff 2 (£l £ ) n L°°(fl £ ). Now setting r = e l , we can choose first k sufficiently large, 
depending on d, n and then / depending on d and k so that F £ is a contraction in the ball B r for e = £j 
sufficiently small. A solution of (| 1 05|) can be then found using the contraction mapping theorem and its 
properties follows from the construction of Uk. e - This concludes the proof. ■ 
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